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■ An exactly soluble non-perturbative model of the pure gauge QCD is derived as a weak coupling limit of the lattice 
theory in plaquette formulation [fij . The model represents QCD as a theory of the weakly interacting field strength 



ON 

oo 



fluxes. The area law behavior of the Wilson loop average is a direct result of this representation: the total flux 
through macroscopic loop is the additive (due to the weakness of the interaction) function of the elementary fluxes. 
The compactness of the gauge group is shown to be the factor which prevents the elementary fluxes contributions 
from cancellation. There is no area law in the non-compact theory. 



! 1 Introduction. 

> ; 

It is well understood that the confinement of quarks can be described by the non-perturbative QCD 
\ only. The discovery of asymptotic freedom || within perturbative approach to the Yang-Mills theory 
unambiguously peaks QCD as a genuine model of strong interaction. However, the growth of the effective 
coupling constant towards infrared region, which is the opposite prediction of asymptotic freedom, makes 
Q\ [ perturbative theory unreliable in infrared. Since we cannot make sense of the path integrals except the 
gaussian case, the derivation of the non-perturbative solution of the Yang- Mills model, which is equivalent 
■ to the practically impossible summation of all Feynman diagrams, cannot be considered as a realistic 
Qh! goal. 



The unique known opportunity to solve the problem non-perturbatively is provided by the Wilson 
lattice model || which is exactly soluble^] for any coupling constant (i.e., non-perturbatively) on any 
finite lattice. In two-dimensions, for example, the lattice model easily gives an exact solution ||] in 
continuum. In D > 2 the problem is much more complicated. However, since the lattice model already 
contains all of the non-perturbative information, the remaining problem of taking continuum limit of the 
lattice results, in spite of being also a non-trivial task, is essentially more attractive. 

It is very important, in this respect, to formulate the lattice model in the most adequate terms. We 
shall argue that such formulation is the plaquette formulation performed for D = 3 in Q (this was given 
in the form different from that of the pioneering works ||). This formulation is the starting point for 
present analysis, and will be briefly reviewed in the next section. The motivation (given also in |IJ) is as 
follows. 

The partition function of the Wilson pure gauge model in D dimensions can be written as 



Z = J UdUiU dU P^ triUp+Utp) K U ^ Ul ) ' w 



where A G is the bare coupling constant, I and p denote links and plaquettes of D-dimensional lattice, Ui 
is the unitary matrix (U(N) or SU(N)) attached to l-th link. The products of C/j's are ordered according 
to the geometrical order of links. 



At least, numerically. 
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The gauge-invariant 5-function is 



5(U,V) = J2xr(tf) X r(V) , (2) 
r 

where r is an irreducible representation with the character Xr(U) and dimension d r = Xr{I)- 

The only non-zero observables in this model are invariant ordered products of link variables along 
closed loops, such as 

w &) = <4 tr n u i) ( 3 ) 

(C is the closed contour). It is understood for a long time that the plaquette variables || [7j are relevant 
to solution of the quark confinement problem ||] . As argued in || , the area-law behavior of (||) , 

W(C) ~ e~ aA , (4) 

where A is the area of the minimal surface S bounded by C, and a is the positive parameter (string 
tension), means confinement of quarks, since it corresponds to the linear potential between colored objects. 

The practical use of the idea of dealing with only the loop variables has been first explored by 
Halpern and Batrouni Q in their field strength and plaquette variable formulation of gauge theory and 
by Makeenko and Migdal Q in relation to the loop equation (see also and references therein). However, 
in spite of all efforts and the progress achieved in these directions, the solution to the confinement problem 
has not been found. 

We follow here a variant (l| of approach || which is based on such formulation of the model in which 
the only independent variables of the model are the plaquette matrices U p . This reformulation requires 
calculation of the integral 

/n^rK^n^) (5) 

J i p lep 

in (||), which is responsible for the interaction between plaquettes. 

Physically, the plaquette matrices are nothing but extended to the compact group manifold field 
strength fluxes (just like link matrices are the group-extended vector-potentials). Confinement, or area 
law, emerges as the consequence of the weakness of the interaction between fluxes - the total flux through 
macroscopic loop is the additive function of the elementary fluxes (plaquettes), i.e. proportional to the 
minimal area. The additiveness of elementary fluxes is the necessary but not sufficient condition for the 
confinement: one has to find a reason why these contributions are accumulated and not canceled in the 
total flux. We will find that the reason is the compactness of the gauge group. We will also show that the 
additiveness is indeed emerges in the A G — > limit (which is the only interesting limit as it corresponds 
to the continuum). In this limit, as it is easy to see from (|l]), the leading contribution comes from the 
saddle point U p = I, and indeed a certain factorization property has to emerge. This leads to a new 
model (described below) where the average (||) takes the form []: 

W ( C ) = (U tJ W L ) (6) 

pes 

where S is arbitrary surface bounded by C. In the model we will present, this average is surface- invariant. 
We emphasize that the average (g) considered in the model (|l|) is not surface- invariant, except abelian 
case (where (0) coincides with (||)). 

The providers of required surface- invariance are the ^-functions of the integral (|5|). In paper this 
integral has been computed exactly (for D = 3), thus, providing the starting point for present analysis. 

The paper organized as follows. In Section 2, we briefly repeat the plaquette formulation of |l]] and 
write the result as a statistical model. In Section 3, we suggest factorization of characters in the weak 

2 Such average, called "filled Wilson loop", first considered in |9| in a different model. 
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coupling limit of the lattice theory and derive the effective model of QCD in arbitrary D. In Section 
4, we discuss compactness and suggest a new method to compute the averages in the compact theory. 
We compute the Wilson loop average and string tension in general form and give a new criterion for 
confinement. We also show that in the mean field approximation our solution indeed gives confinement 
and reproduces the known results in both strong coupling and perturbative regimes. In Section 5, we 
summarize and discuss the results. We claim that our results mean the rigorous proof of confinement in 
QCD. In Appendix, we give some useful group-theoretical formulas relevant to the heat kernel analysis. 

2 Lattice QCD3 as a statistical model. 

After taking integral (||) and substituting it into (|l|) one obtains JjJ] the model which is defined on the 
two-dimensional infinite genus lattice formed by hexagons h as shown in Figjlj. 




(a) (b) 

Figure 1: (a) Thin lines correspond to the links of the original 3d lattice. Thick and dotted lines are links of the resulting 
2d lattice. This 8-cube fragment is described by function F in (^) . Integration goes over dotted links only, (b) The fragment 
of the resulting (smoothed) surface, corresponding to 32 cubes of original lattice. 

Here, we again temporarily introduce link unitary matrices, which are the plaquette matrices of the 
original model ([l]). The partition function is 

Z=fUdU l e£*C«W/)Jl5( J] U h , J] U h ) . (7) 

1 h heh heh 

The <5-function is again given by (|2|). It is important that the order of hexagon's links in the product 
(argument of (5-function) is fixed (as hfahfahfa) so that each character contains either l\ or fa-type link 
matrices (in Fig. 1(a), they are denoted by either thick or dotted lines respectively). A fragment of the 
resulting 2d lattice is shown in Fig. 1 (b) . The whole lattice is easy to imagine as obtained from the regular 
3d lattice (the lattice spacing is doubled) by replacing links by the tubes (handles) and vertices by the 
smooth connections of tubes. 

Besides this two-dimensional property of the model, there is another important consequence of the 
plaquette formulation which has not been emphasized in paper 

Namely, we notice that the partition function (0) can be written in the form of the statistical model 
of the integer valued A-component field r: 

Z = Y,U F ^o;r 1 (j),...,rsU)), (8) 

{r} j 

where {r} means a sum over all configurations of r's with each r being assignment of hexagon; index j 
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labels the 8- hexagon configurations shown in Fig. 1(a) and described by 

12 8 

F(X ;n,...,r 8 ) = [jjdUk e £ tr ^ +£/ *> f[ Xn (U itl U it2 U i>3 ) . (9) 

J h=l i=l 

The integral is taken over 12 links corresponding to the dotted lines in Fig. 1(a). Only these links enter the 
characters under integral. Index i in the second product of (||) labels the hexagons in the configuration. 
One should emphasize that the surface of Fig. 1 (b) separates the three-dimensional space into two identical 
subspaces, and j labels the configurations from both subspaces (otherwise there would be no interaction 
in the model). 

This representation unambiguously demonstrates the local character of the interaction between the 
plaquette matrices U p (we remind that the matrices of @ are originally the plaquette matrices U p ). The 
integrals over them enter the partition function via the function F only. Nothing similar to this result 
can be observed in the link variables formulation of the model. This locality means that the interaction 
between plaquette matrices U p gives rather weak influence on a long-distance (continuum) phenomena 
which are, instead, defined by the interaction between iV-component fields r (irreducible representations). 
The contribution due to U p s affects only the form of the function F, or eventually, the form of the action, 
and has to be irrelevant in the continuum, especially for such a phenomenon as confinement. 



3 Weak coupling limit of the lattice QCD. 

It is well known that the continuum limit corresponds to the weak coupling, A Q — > 0, limit of the lattice 
model. This has to appear as a result of the renormalization-group analysis on a lattice, but also can be 
used in advance. Namely, we take the weak coupling limit at the fixed lattice size, while the refinement 
of the lattice will be postponed for the later stages. 

In the A D — > limit, the leading contribution to the function F comes from the abelian saddle point, 
Up = 1. In the vicinity of this point, the model still preserves its non-abelian nature since the restriction 
Up = I for the plaquette matrix (which is the product of link matrices, U p = Y\i£ P Ui), imposes a slight 
constraint on each link matrix Ui, and leaves its non-abelian property unbroken. 

Technically, near this saddle point, the character of the product can be replaced by the product of 
characters, 

Xr(ABC) -» d- 2 Xr (A)xr(B)xr(C) , (10) 

where factor d~ 2 provides a proper normalization. Several initial terms of the saddle-point expansion 
of F are exactly given by the substitution fllQ|). In abelian theory, in particular, both sides of (10) are 
identical to each other. 

After replacement (|H]), D = 3 partition function takes the form 

2=En<n/p. en) 

{r} c p 

where sum goes over all configurations of r's, product goes over all cubes c of the original lattice, and 
interaction between neighboring cubes c\ and c%, on their common plaquette p is defined by 

f p D = 3 = fdU ^(U^XrAUjXr.m (12) 

Here and below, we do not distinguish in notations between cubes and representations attached to them. 
For example, n means representation attached to the cube c\. 

Applying ([ji]) to the Wilson loop average (||), we realize that it takes the form of Such an 
average, if considered in the model (jl|), is not surface-invariant (except the U(l) case). This is easy to 
see in the large- N limit where this average factorizes into the one-plaquette averages W a with the result 
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W{C) = (S is arbitrary surface). The reason for this is simple: factorization (|To| ) deforms (except 

the U(l) case) the (^-functions which are the only providers of the surface-invariance in (|l|). Nevertheless, 
in the model ( [TT] ) (which is obtained from {JlJ) by the same deformation (|T0|)), as we will check below, the 
average (P) is surface-invariant. 



The physics behind the replacement (10) is quite simple. Let us remind that the plaquette variables 



are the group-extended field strength fluxes. The replacement (10) manifests the abelian character of 
interaction between them. Let us emphasize that this property is clearly presents already in the original 
model (H). In the weak coupling limit, this property becomes more transparent. It is this property which 
provides the first important ingredient of the area law, as it leads to the additiveness (in terms of the 
algebra) of the elementary fluxes. To make this qualitative picture more precise, we will perform the group 
Fourier-transform of the model, i.e., complete the formulation in terms of the irreducible representations 
r only (r is the Fourier-image of the unitary matrix U). 

In arbitrary D, even without a step- by-step repetition of the procedure of [||, it is easy to understand 
that after the factorization (10), the partition function takes the same form (|TT|) , with the only difference 



being in the function f p . In D = 4, f p is the integral of four characters (instead of two in D = 3) 
corresponding to four cubes c p sharing plaquette p. In the arbitrary- D case, f p is the integral of 2D — 4 
characters corresponding to 2D — 4 cubes sharing each plaquette: 

= JdU e£ tr(C/+C7t) 'ff^f^ • (13) 

C v = l Tc P 



As a remark, it is fruitful to understand how the model (|Tl|) emerges in the link-variables formulation. 
Considering cube as a frame (set of its links), one writes its two-dimensional functional of boundaries 
according to the general formula pR: 



m 



Jd^e-W)!!^), (14) 

r p=l 

where Ci{r) is the quadratic Casimir eigenvalue. For the cubic frame, m = 6, and .A = 0. Thus, ( |i~4|) 
takes the form 

E d r 4 UXr(U p ). (15) 

r p=l 

Then, the product of these factors over all cubes, integrated over independent variables U p 's with the 
weights exp j-tr (U p + Up~), i.e|], 

Z=[ \{dU p e^+^nEC A Xr„(EW , (16) 

p c r c p c =l 

exactly coincides with the partition function (pd|). 

Now we compute the function f p . Irreducible representation r is parametrized by the Young table 
parameters (representation's highest weight components) n M (/i = 1, N) with the dominance condition 
n\ > n<2 > ... > n N . We switch to ha = — \i + ./V (h\ > hi > ... > h N ). The formula for dimension of 
representation is 

d r = ^§-y A(h) = H(h,-K) (17) 

(A(/i) is the Vandermonde determinant, h ( j 1 = N — fx), and the character is: 

c" 

A(e**) " r A(e 



»(£0-%^-+^St/^^», (18) 



3 This model, in a different context, was considered in [hoi 
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where e^ are eigenvalues of U. The second equality in (|T^) is due to the Itzykson-Zuber formula [11]. 
Substituting this formulas to the equation (12), we have 



r 



D=3 



det M „ I h/j (ri ) _ hv (r2 ) (2N/ X ) 



where I n is the modified Bessel function. Replacing I n {2N / \ Q ) by its A c —* asymptotics, 



we write (19) as 



D=3 



„(_) „ e" - 



A(hx)A(h 2 ) 



(19) 



(20) 



(21) 



In the Appendix, we show that substitution ( p0[ ) is equivalent to the replacing Wilson action exp ^-tr (U+ 
W) by the heat kernel f§: J2 R d R 

point (A — ► 0) expansion around (j) = (U = I) in(0). 

In arbitrary D, substituting (|i~7|) and (|l8|) to (|13|), we have 



e 2 n C2 ^Xr(U). This replacement is also equivalent to the saddle 



r 



A' 



n(< 

!/=l 



^ e ^ cos< MA 2 (0) 



A(0) 



A( e' 1 



2(D-3) 



2D-4 

n 



itr 4>V Cp h Cp Vj p 



(22) 



The factor A( &i 
expansion around 



AW 



(emerges in D > 3) gives no contribution to the leading order of the saddle-point 
= 0. Ignoring this factor, we take the gaussian integral over cj> and obtain 



. 2D-4 



fp=J Y[ dV Cp e~^ tr{ ^ VcphcpVctp)2 



(23) 



In D = 3, this reproduces (pTT) . In D 



ft 4 = n 



e 2N lrn j 



= 1 

4, this gives 
6 



/ n ^e> w «n^ n ^) 

(y)=i fc=i (w)efe 



(24) 



{ij)ek 



(we have switched to the notation hj = h(cj)). Geometrically, formula ( f24[ ) can be represented as a 
tetrahedron with vertices j, links (ij) and triangles k, with corresponding assignment of matrices hj, 
Uij = ViVj and the conditions at the triangles. 

As one can see from this example, the reduction of ( |23| ) to the eigenvalues, like (0), is highly non- 
trivial in D > 3. Exactly the same problem appears in the matrix models of 2d quantum gravity in 
the case when embedding target space corresponds to the physical situation of central charge c > 1. 
In spite of the efforts made in recent years the solution to such problems is unknown, though some 
important physical information can be extracted in the simplest physical case of only one constraint (this 
corresponds to c = 1), considered in [|12] (see also revie w [|l~3| ] and references therein). 

In spite of the highly non-gaussian character of f p (g^) in terms of eigenvalues h, as one can see in 
examples of ( ^l| ) and ( |24| ) , f p is the gaussian function in terms of the full hermitian matrices H = VhV^ . 

We shall use f p in the general form (23). The resulting partition function is 



{h} c p c p 



(25) 



(cp is the cube containing plaquette p). The compactness of the gauge group is reflected in the integerness, 
or rather discreteness, of eigenvalues h. 
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4 Compactness, loop averages and string tension. 



As we already noticed, the Wilson loop average @ takes the form (g) in the model (p5|). At the plaquettes 
p £ S, f p has to be replaced by 



The eigenvalues of matrix -ff(f) are the fundamental representation components; tr H 2 {i) = C*2(f) where 
C*2(f) is quadratic Casimir of the fundamental representation: C^f) = -/V for U{N) and C*2(f) = A — 
for SU(N). Then, 

= (II f ) = e-^^(II II e-^ H ^)) , (27) 
pes pes c p 



where averaging (...) is understood in the sense of (p5|), A is the lattice area of S. 

The average (|27| ) is manifestly surface- invariant. To check this, one chooses another surface S' and 
makes the shift H —> H — H(f) in the 3- volume connecting surfaces S and S' (in D = 3 this is simply 
a volume bounded by the compact surface S U S'). The result is again ( p7| ) but with the surface S'. 
The global minimum of the action corresponds to the minimal surface. However, one cannot vary over 
integer- valued variables. 

We apply now the generalization of the idea [14] realized in two-dimensional large- A theory. In the 
non-compact theory, h is continuous variable. Then, the sum over /i-configurations in ( |25| ) turns into the 
path integral / T>h(x) with the measure T>h = A 2 (h) JJ k dh^. The dominance constraint is automatically 
satisfied due to the antisymmetry of the Vandermonde determinant. In the compact theory, we still can 
redefine the nature of variables assuming that the effect of compactness can be accumulated into the new 
coupling constant A. The definition of A is 



(28) 



{h} c p 
In other words, the model takes the form 



z = J n m ( c ) n / n e ~ 2> ^ Ecp Ucv > ■ ^ 



Since in the theory defined by (^9|) one can vary over H, the minimal area (which corresponds to the 
maximum of the integrand) becomes distinguished. We have 

W{C) = e-^^-^ JJ jj e -^ to ^^( f )) non _ compact = e~ L ^ 2m ^W'(C) , (30) 

where W'(C) is the average in the non-compact theory with coupling constant A defined by (^8|). 
In D = 2 case of the infinite plane, 

W'(C)= e~^ A => W{C)= e — dh A (31) 

for any N, including abelian theory, in agreement with the known result [||. It is also not difficult to 
check that in the topologically non-trivial D = 2 cases, one can reproduce the results of 0. 

In D > 2, W'(C) contains no area term. According to the definition W(C) ~ e — mln , equation ( |30| ) 
gives the dimensionless (since we have a dimensionless area) string tension 

2 N~ ' (32) 
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where A is defined by equation (p8|). In general, A is the function of all parameters of the theory: A Q , N 
and D. 

As an example, we demonstrate that in the case of the compact QED3, (^) reproduces the result of 



[15|. Equation (|28|) takes the form 



E e -¥E w K-^) 2 = f Udhc e -tE< w (>^) 2 (33) 



{n} 

In the mean-field approximation we have 



f dh e"^ 2 = Y, e-^" 2 = fdh e~^ h2 ]T e"^™ 2 (34) 



(the second equality is the Poisson resummation formula; m is integer). Thus, 

1 1 



and 



E e "^ m (35) 



in 



°vm = ^7T^ = TT f 1 ^ 2 : ) , (36) 



2 (Em e"^ m2 ) 2 



i.e., gives both small- A result of [15]: c ~ A e~^, and the known large-A D behavior: a ~ 4f. 

It is also clear that the qualitatively similar A D —* and A D — > 00 behaviors take place in the 
non-abelian case and thus coincide with the known results of perturbative analysis and strong coupling 
expansion. 



Similarly to (30), the loop average in representation R is 

Wr {C) = ( Xfl( ^ (C)) ) ~ e- a Wr A ^ . (37) 
This gives a possibility to compute the arbitrary average, 

^) = m = ew^). (38) 

R 

where 

U = J dUt(U) X n(tf) (39) 

is the Fourier coefficient of £,(U). 



5 Summary and discussion. 

Thus, taking the A Q — > limit of QCD in plaquette formulation Q, we derived the non-perturbative 
model (p5|), which is formulated in terms of observables only (loop variables) and therefore is much more 
suitable than the model ([!]) for the analysis of physics of QCD. This model is shown to be exactly soluble 
in some particularly interesting cases relevant to the continuum. Although the model (|25|) is obtained as 
the A — > limit of the Wilson model (|l|), we suggest that it also describes the lattice QCD in the non- 
empty interval A Q £ [0, A*]. By definition, A* is the upper limit where the model ( p5| ) is still equivalent 
to ©. In the 17(1) case, for example, A* = 00 since (|25| ) coincides with (|l]). 

We computed continuum Wilson loop average and obtained the area law. Since equation (^8[) gives 
A D > A, the string tension (^) is non- vanishing. We consider this result as eventual proof of confinement 
in QCD. 
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Qualitatively, the picture of confinement can be given in terms of elementary fluxes, which are the 
plaquette variables. Then, QCD is a theory of weakly interacting fluxes. The total flux through any 
surface is the additive function of elementary fluxes, which provides the necessary ingredient for the area 
law. The interaction between fluxes provides the surface-invariance of the averages like @, which is the 
eventual result of the global gauge invariance. 

The reason which turns this additiveness into the area law (and thus provides the sufficient ingredient) 



is compactness of the gauge group (reflected in A G > A solution of the equation (28) and thus into the 
non- vanishing string tension). A solution of the equation is non-trivial problem in general. Here, 
we only demonstrated that it gives proper behavior for the string tension both in the A D —* and in the 
A — > oo limits, and thus provides A G > A for any A G . 

The additiveness of fluxes takes place also in the model (|l|) and becomes more transparent in the model 
(25) due to the factorization of characters (10). The factorization also provides the crucial simplification 



of the model (as one might expect in the continuum limit) and its solubility. As we already mentioned, 
the reason why the model remains non-trivial under this replacement lays in the nature of plaquette 
formulation: the condition U p = J\i Ui = I imposes a slight restriction on matrices XJ\. Thus one learns 
that even within the gaussian model the non-trivial results can be obtained. Recently, similar idea was 



argued by I.Kogan and A.Kovner [16]. 
There are several aside remarks. 

There is a complete algebraic similarity with the principal chiral field model (PCF) defined in the lower 
dimensions than QCD. The above considerations can be easily carried on for PCF: the PCF/) partition 



function is given by QCD^-i (25) where the cubes and plaquettes should be respectively replaced by 
the plaquettes and links []. The Wilson loop (|3|) corresponds to the two-point PCF correlator G xy . The 
plaquette formulation pj becomes the link formulation for PCF, and G xy takes the form of (p7|), where 
plaquettes of the surface S should be replaced by the links of arbitrary curve C xy connecting points x 
and y. To prove that another curve C' xy in D = 2 produces the same result, one makes the constant 
shift H — ► H + H(i) in all plaquettes of the surface enclosed by C xy U C xy . To prove the invariance in 
D = 3, one has to make such a shift in all plaquettes of the arbitrary surface enclosed by C xy U C'. We 
actually used the analogy between QCD and PCF to check the surface-invariance of the loop average 



in QCD 4 . 

It is easy to recognize in D = 3 version of (|2^) the Kazakov-Migdal model |17]]. Even though the 
coincidence takes place only in D = 3, the observation is quite peculiar. The idea of |l7[] to interpret 
the diagonalizing unitary matrices U xy as a prototype of the QCD gauge field was dismissed due to 
the observation [||, |l^] that [/^-built averages are trivial due to the local Zjq symmetry of the model. 
However, as we now understand, the KM master field (i.e., solution of the large- N saddle point equation 
for the partition function) is indeed true QCD master field in D = 3 (and only in D = 3), while U xy has 
nothing to do with the QCD gauge field (the actual gauge field of QCD is represented by its Fourier- 
images, h x fields). Thus, the KM master field could be directly used in the large- N limit of the model 
(29) in D = 3. Unfortunately, this is not found]^. 

It is still a challenging problem to extend the present analysis to the full version of the model which 
includes the matter fields. 

The highly non-trivial problem is the renormalization-group analysis and computation of the non- 
perturbative beta-function. I suppose this can be achieved in the computational direction provided by 
the equation (pq). 

Acknowledgements. I am grateful to D.Boulatov, M.Halpern, I.Klebanov, I.Kogan, A.Kovner, 
M.Peskin, A.Polyakov, M.Shifman and W.Taylor for fruitful discussions. 



4 We imply only algebraic similarity. Regarding some dynamical properties, this might be different. For example, QCD4 
in some aspects is rather similar to PCF2 than to PCF3. 

5 The solution claimed in jl9| disagrees with that of E3|, and the discrepancy has not been resolved. 
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Appendix: Heat kernel and A G -expansion. 

In this Appendix we give some useful formulas and establish a connection between the technique of 
Section 3 and the heat kernel technique. 

In derivation of (pT[|), the Bessel function has been replaced by its A Q — > asymptotics. We now 
demonstrate that this is an equivalent of the replacement of the Wilson action exp £tr(c7 + [/t) by the 
heat kernel 0, 

J2d R e-^ C ^Xn(U). (40) 

R 

Expression ( |T2| ) after substitution ( f40| ) takes the form 

f P = -f±rJ:d R e-^ {R) D« r2 , (41) 
where D^ ir2 is the multiplicity of representation R in the tensor product r\ <g> fj, 

= J dU X R(U ] )xrAU)Xr 2 {U ] ) , (42) 
Xri (U)Xr 2 (tf) = J2Kr 2 Xn(U) . (43) 

R 
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The equivalence between (21) and (Ell) can be easily established in their A -expansions. Expansion of 



Itzykson-Zuber integral entering fl2l|) is 

^f,T = fdU e~^ AUBUt = Y tMl f dU tT n AUBU (44) 
A(o)A(6) J ^ n\N n J y ' 

(where a and b are the diagonal matrices of the eigenvalues of A and B). The low orders of this expansion 
can be derived using explicit expressions for powers of traces tr n AUBU via the characters Xr(AUBW) 
and formula 

J dUxMUBU^) = ±- Xr (A) Xr (B) . (45) 



We have 



dct M „ e~— a » 0v -. \o 



1 - i&tr AtrB+ 



(46) 



A(a)A(6) ~ 1 N 

( Ntv 2Atv 2 B + Ntv AHv B 2 - tr Mtr B 2 - tr ^ 2 tr 2 #) + 0(A*J) , 
where to derive 0(X 2 ) term we have used 

X20...oOO = ^(trV + try 2 ), XU0..AV) = ^(trV-tr^ 2 ) , (47) 



and thus, 

,2 



tr 2 ^ = X20...o(^)+Xiio...o(^) • (48) 



The higher order terms can be derived similarly. Expansion (|46D is then to be substituted into (|2l|). 
Expansion of (41) can be derived by technique of |2l]] which is based on the replacement of —C 2 {R) 



by the group Laplace operator trdy, 



f P = e™ tld hrAU)Xr 2 (U^\ u=i . (49) 



Applying tr dfj to both sides of ( |43[) one derives the identity 

-^L-J2 D nr 2 dnC 2 (R) = C 2 (n) + C 2 (r 2 ) - ^Ci(n)Ci(r 2 ) , (50) 

which gives the order 0(A ) of f p . C\{r) is the first (linear) Casimir eigenvalue. To obtain the next order 
we apply tr dfj twice to ([I||) and get 



(51) 



d^r 2 T, R D« r2 d R C$(R) = (C 2 ( n ) + C 2 (r 2 ) - £Ci(n)Ci(r 2 )) 2 + 

^(c 2 (n) - j?cKn)) (C 2 (r 2 ) - ±C 2 {r 2 )) , 

etc. Thus, for (|4l|), one has 

f P = 1 - &(Oj(ri) + C 2 (r 2 ) - f Ci(n)Ci(r 2 )) + 

A(0,(n) + C 2 (r 2 ) - 2 C 1 (r 1 )C 1 (r 2 )f+ (52) 



N ^-i) ( C 2(n) - jjC 2 {n))(c 2 {r 2 ) - jfC 2 (r 2 ))+0(\l) . 
Taking into account explicit expressions for Casimir eigenvalues, 



N 

C 1 (r) = Y / h k =trh, (53) 
k=l 
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C 2 (r) =tvh 2 -R , (54) 

(R = N ^ N 12 ^ and definition of is changed to = n^ — /x + ^^i) one checks order- by-order agreement 
between @ (with the substitution of @) and (||). 

Thus, the heat-kernel Young table parameters play the role of the eigenvalues of the hermitian 
matrix of Section 3, and the subsequent equivalence between heat-kernel and matrix-model approaches 
one establishes by identification: C2(r) with tr/i 2 , C\{r) with tr/i, d r with A(h) etc. The A D -expansion 
(|46D, or (f)2j), shows that / p is non-singular symmetric function of h. 

There are some changes for SU(N) group with respect to the general unitary group U(N). In 
particular, Casimirs of SU(N) can be obtained from Casimirs of U (N) by the replacement — > b! = 

- ^tr h. Thus, for SU{N), 

Ci(r) = 0, (55) 

C 2 (r) =trh 2 - -^tr 2 h (56) 

(we ignore the additive constant R). Thus, hermitian matrix h of Section 3 is the traceless matrix if the 
gauge group is SU(N). 

Correspondingly, for SU(N), the formula (|5T|) takes the form 



^L-J2 D ^r 3 dnCi(R) = (c 2 (r 1 ) + C 2 (r 2 ))\ J ^— [ C 2 (r 1 )C 2 (r 2 ) , (57) 
while instead of (|5^) we have 

I, = 1 - ^(ftW + C 2 (r 2) ) + J,( C2 (n) + C 2 (, 2 )) 2 + ^ff^ 1 ' + £)(*.•) . (58) 



In applications, the generalization of the formula ( p0[ ) to the case of arbitrary number of characters can 
be useful: 

J2d R C 2 (R) f ^n¥^ = E^W + ^E C iWCi(r;). (59) 

R j r 3 j i<j 
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